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ABSTRACT 

The eigenfrequency problem of fundamental vibration mode in a highly inho- 
mogeneous star modeled by self-gravitating mass of viscous liquid with singular 
density at the center is considered in juxtaposition with that for Kelvin funda- 
mental mode in the liquid star model with uniform density. Particular attention 
is given to the difference between spectral equations for the frequency and lifetime 
of /-mode in the singular and homogeneous star models. The newly obtained re- 
sults are discussed in the context of theoretical asteroseismology of white dwarfs. 



1. Introduction 

The term fundamental mode has been introduced in the theory of stellar pulsations 
by Cowling (1941) from hydrodynamical model of a heavy mass of inviscid incompressible 
liquid of uniform density p undergoing free oscillations with nodeless irrotational velocity of 
fluctuating flow. The extended discussion of the eigenfrequency problem of this model, whose 
solution is due to Kelvin (1863), can be found elsewhere (e.g. Lamb 1945, Chandrasekhar 
1961, Bastrukov 1996). An outstanding importance of this fiducial model for theoretical 
asteroseismology is that it sets the standard for analytic study of non-radial pulsations of 
the main-sequence stars (Cox 1980, Unno et al 1989) and serves as an example, according to 
Chandrasekhar (1961), of "at least one problem for which analytic solution can be found" and 
which illustrates " the type of difficulties one must confront in the other problem" , because 
" in most instances, the problems become of such complexity and involve so many parameters 
that elementary methods of solution seem impracticable". 
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The purpose of this work is to explore some peculiarities of the fundamental vibration 
mode in the inhomogeneous stars. In approaching the topic like this it seems better to start 
with a brief outline of assumptions lying at the base of the /-mode frequency computation in 
the homogeneous liquid star model in which the uniform equilibrium density is not altered, 
5p = 0. The obvious consequence of this assumption steaming from the continuity equation 

Sp = _V(p5v) = -(<Jv V)p-pV-5v (1) 

is that the velocity of flow oscillating in a spherical mass of incompressible liquid of uniform 
density is described by the potential vector field of the form 

V-5v(r,t) = 0, <fv(r,t)=Vx(r,*), (2) 
V 2 X (r,t) = 0, X (r,t) = A e r e P e (cos6)a(t). (3) 

The last equation exhibits nodeless character of the velocity flow as a function of distance 
from center to the surface of the star. It is this feature of oscillating flow is regarded as the 
major kinematic signature of fundamental vibration mode. 

In this paper we relax the basic assumption of the homogeneous liquid star model about 
uniform equilibrium density, preserving the above major kinematic signature of /-mode. 
Specifically, we consider admittedly idealized a highly inhomogeneous liquid star model in 
which the non-uniform density profile has singularity at the star center of the form 

P^ = Ip°(j) 1/2 - ( 4 ) 

One of the most conspicuous features of this density profile is that the total mass of such a 
star 




is finite and identical to that for the total mass of the canonical homogenous liquid star model. 
In somewhat different context the model of inhomogeneous liquid star with similar density 
profile has been briefly discussed by Clayton (1986). This curious feature of the model, which 
from now on is referred to as the singular star model, is interesting in its own right because 
it permits too analytically tractable solution of the eigenfrequency problem for /-mode. 
Perhaps the most striking feature distinguishing vibrational behavior of inhomogeneous from 
homogeneous one is that the nodeless non-rotational vibrations of a spherical liquid mass 
of non-uniform density are of substantially compressional character, as it follows from the 
continuity equation. Therefore, one of the prime purposes of our study is to find out how this 



- 3- 



striking distinction between the density profiles of homogeneous and inhomogeneous models 
is reflected in the frequency and lifetime spectra of fundamental vibration mode. 

The paper is organized as follows. In section 2, a general mathematical treatment of 
the inhomogeneous star vibrations in the fundamental mode on the basis of the Rayleigh's 
energy variational principle is outlined. In section 3, the detailed analytic computation of 
the spectral equations for the frequency and lifetime of /-mode in the singular liquid star 
model is presented followed by their comparison with those for the Kelvin fundamental mode 
in homogeneous star model which is the subject of Section 4. The newly obtained results 
are highlighted in section 5 and briefly discussed in the context of current development of 
theoretical asteroseismology of white dwarfs. 



2. General equations of nodeless stellar vibrations in fundamental mode 

The state of motion of flowing stellar matter under the action of forces of buoyancy, 
gravity and viscous stress is uniquely described in terms of five dynamical variables, to wit, 
the density p(r,t), three components of the flow velocity v(r, t) and the pressure p{r,t) 
obeying the coupled equations of fluid-mechanics and Newtonian universal gravity 

dp dvk d d d .g. 

dt ^ dxk ' dt dt Vk dxk ' 



d 
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dt 


= m +Vk 


dU 


, dir ik 



^dvj = dp | | „.. tK 
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dp T±pdp dv k dhip 

dt p dt dxk amp 

V 2 [7 = -AnGp. (9) 

The equation for density p is the condition of continuity of flowing stellar matter. The Navier- 
Stokes equation for the velocity flow expresses the second law of Newtonian dynamics for 
viscous liquid and equation for pressure is the condition of adiabatic behavior of stellar 
matter; Ti stands for the adiabatic coefficient of gaseous pressure in the star. This latter 
condition means that in the process of motion, like vibrations, the time scale of energy 
exchange between infinitesimally close elementary volume of liquid is much longer than 
characteristic period of oscillations (e.g. Huang 1987, Kawaler & Hansen 1994). The tensor 
of Newtonian viscous stresses is given by 



2 \ 1 

2 



n ik = 2r] v ik + ( C - o V I v n $ik, v ik = 7T (VfcVj + ViV k ) ■ (10) 
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In computing frequency of fundamental vibration mode, the equilibrium density profile p, 
the pressure P c at the star center as well as the transport coefficients of stellar matter, the 
shear rj and the bulk ( viscosities, are regarded as input, in advance given, parameters. 

The potential of self-gravity U(r) and pressure p(r) in motionless, v = 0, state of 
hydrostatic equilibrium, are the solutions of coupled equations 

Vip(r) = p(r)Vi U(r), V 2 U(r) = -AnGp(r). (11) 

The gravity potential U is determined, in effect, by Poisson equation for internal potential C/j 
and Laplace equation for external one U e supplemented by the standard boundary conditions 
of the continuity of these potentials and their normal derivatives on the star surface 

V 2 U t (r) = -AnGp(r), r < R, V 2 U e (r) =0, r > R, (12) 

U t (r) = U e (r) , djm = djUr± (13) 
r=R ar ar r=B, 

The general solution of equation for pressure is specified by standard boundary condition of 
stress-free surface p(r = R) = 0. 

The equations of linear oscillations of stellar matter about stationary state of hydro- 
static equilibrium of a star with non- uniform equilibrium density, p = p(r), are obtained by 
applying to ©-© the standard procedure of linearization 

p — >• p{r) + 5p(r,t) p — > p(r) + 5p(r,t) (14) 
Vi^Vi + 5vi(T,t), [vi = 0], U ^U(r) + SU(r,t). (15) 

As was stated, in this work we focus on the regime of irrotational vibrations in which the 
velocity field of fluctuating flow subjects to 

V ■ 5v = 0, V x 5v = 0. (16) 

Given this, the continuity equation takes form 

5p = -(6v k V k )p. (17) 

The liner fluctuations of the flow velocity are governed by linearized Navier-Stokes equation 

pSvi = -Vi5p + pVi5U + SpViU + V k 5-K ik , (18) 

5-K ik = 2r)5v ik , 5v ik = ^(V^ + V^). (19) 
From (|S]) it follows that rate of change in the pressure is controlled by equation 

Sp = -(v k V k )p. (20) 
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Fluctuations in potential of self-gravity inside the star 5U = 5Ui, caused by fluctuations in 
density Sp, subject to the Poisson equation 

V 2 5U = -AnG5p. (21) 

The energy balance in the process of oscillations is controlled by equation 

Wt S dV = ~ I ^ 5vkVk)5p - P( 6v ^^ 6U - Sp(6v k V k )U - 6n ik 5v ik ] dV (22) 



which is obtained after scalar multiplication of ffl8l) by 5vi and integration over the star 
volume. To compute the eigenfrequency of vibrations we take advantage of the Rayleigh's 
energy variational principle at the base of which lies the following separable representation 
of fluctuating variables 

Svi(r,t) = ai(T)a(t), 6U(r,t) = <p(r)a(t), (23) 

Sp(r, t) = p{r) a(t), p(r) = -(a k (r)V k ) p(r), (24) 

5p{r,t)=p{r)a{t), p(r) = -(a*(r)V*) p(r), (25) 

5ir ik (r,t) = 2r]5v ik = 2r)a ik (r)a(t), (26) 

Sv ik (r,t) = a ik (r)a(t), a ik (r) = -(Vja fc (r), +V fc aj(r)). (27) 

Hereafter ai(r) stands for the time- independent field of instantaneous material displacements 
and a(t) for the temporal amplitude of oscillations. The key idea of such representation is 
that it transforms equation of energy balance (f22j) into equation for a having well-familiar 
form of equation of damped oscillations 

dS Ma 2 ICa 2 ^ Va 2 

_ = -2* £= — + —, T=— , (28) 
Ma + Va + JCa = 0, (29) 

M = J p(r)a k a k dV, V = 2 / r](r)a ik a tk dV (30) 

K = J [{a k V k p{r)) (a k V k U(r)) - p(r)(a k V fc )0(r) - (a k V k ) (a k \7 k )p(r)}dV. (31) 

Here £ is the energy of free, non-dissipative, oscillations and T is the dissipative function of 
Rayleigh describing their damping by shear viscosity of stellar matter. The solution of ( 1291) 
is given by 

a(t) = a exp(— t/r) cosu;(r)t, (32) 
u\r)=u 2 [l-{ur)-\ u 2 = ^ r = ^ (33) 



- 6- 



where oj{t) is the frequency of dissipative oscillations damped by viscosity, ui is the frequency 
of free oscillations and r is their lifetime. Thus, to compute the frequency and lifetime one 
need to specify all variables entering integral parameters of inertia M., stiffness /C and viscous 
friction T>. 

We start with the field of instantaneous displacements a k (r) which is the key kinematic 
characteristics of /-mode. In the star undergoing irrotational oscillations in the /-mode, 
the shape of an arbitrary spherical surface takes the form of harmonic spheroids which are 
described by 

r(t)=r[l + a e (t)P e (()\, ( = co S (34) 

where Pe(() is the Legendre polynomial of the multipole order £ specifying the overtone 
number in fundamental vibration mode. In the system with fixed polar axis the potential 
field of velocity is found from Laplace equation supplemented by boundary condition that 
radial component of velocity on the star surface equals to the rate of the surface distortions 



taking the shape of harmonic spheroids 

V k 6v k = 0, 5v k = V kX , V 2 x = -> x = A E r e P e (()a(t), (35) 

Sv r \ r=R = R(t), R(t) = R[l + ae(t)P e (Q]. (36) 

Taking into account that 5v(r,t) = a(r)d(t) one has 

a(r) = ^V/P,(C), ^ = ^=2- (37) 
In terms of this field, the linearized equations for the density 5p and for the pressure 5p reads 

5p(v,t)=p(r)a(t), p(r) = -(a • V) p(r) = -^(CK^ 1 ^, (38) 

5p(r,t)=p(r)a(t), p(r) = -(a- V)p(r) = -^(CK^ 1 ^ (39) 



where p(r) and p(r) are the density and pressure of gravitationally equilibrium, hydrostatic, 
configuration. 

To compute variations in the potential of gravity one must consider two equations for 
internal SUi and external 5Ui potentials 

V 2 <JC^(r, t) = -4tt G5p(r, t), V 2 0,(r) = -4tt Gp(r) (40) 
V 2 5U e (r,t) = 0, V 2 e (r) = O. (41) 

The outlined energy variational method provides a general framework for computing the 
frequency of /-mode in a inhomogeneous Newtonian liquid star with arbitrary form of non- 
uniform density profile. 
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3. The singular star model 



As was stated, we use the term singular star for a self-gravitating mass of viscous liquid 
with the non-uniform and singular in the star center density profile given by equation 
whose total mass is identical to that for homogenous star model. For such a singular star, 
the equilibrium potentials and fields of universal gravity are given by 

1/2 



V 2 f/,(r) = 
Ui(r < R) 
gi(r < R) 



-AnGp(r), 
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(42) 
(43) 



g e (r >R) = -VU e 



GM 



-r. 



The hydrostatic pressure obeying the boundary condition of free surface p(r) 



p(r) 



lOvr 



Gp 2 # 



P, 1 
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Pr 



IOtt 



5 GM 2 



(44) 
= reads 

(45) 



It is remarkable that the pressure at the star center, where the density has singularity, is finite 
and its radial profile is linear function of distance from the center to the star surface. The 
pressure in the center is defined by equation of state of the stellar mattero With in advance 
given equation of state for P c the rightmost identity in (j4"5"l) is considered as definition of the 
star radius. The internal gravitational energy is 

11 GM 2 . _„ _ 3GM 2 

(46) 



2 / PsUidV 



1.02W, 



5 R 



18 R 

where W is the total gravitational energy of homogeneous star of equivalent mass M. To take 
into account the compression effect of self-gravity on mechanical property of star matter we 
adopt the suggestion of work (Bastrukov et al 2007a) that radial profile of viscosity coefficient 
is identical to that for the equilibrium pressure profile, that is, of the form 

r 

R, 

where u c is the shear viscosity in the star center which along with p s and P c are regarded as 
input parameters of the model. 



Hr) =u e (l- 



(47) 



1 In white dwarfs the central pressure is identified with the pressure of degenerate Fermi-gas of ultra 
relativistic electrons: P c = Pf = Kn'J 3 . The remarkable feature of the emergence of these end products of 
stellar evolution (e.g. Hansen & Kawaler 1994) is that the density at the center of a white dwarf tends to 
infinity as the mass of collapsing star approaches the Chandrasekhar limiting mass (e.g. Phillips 1994). 
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3.1. Exact solution of Poisson equation for variations of self-gravity potential 

in singular star 

Having denned the equilibrium profiles of density p(r), the pressure p(r), the shear 
viscosity profile z/(r) and knowing the field of instantaneous displacements a(r) we are able 
to compute the inertia Ai and viscous friction V. However, in order to compute the stiffness 
/C we must calculate fluctuations in the potential of self-gravity 5U(r,t) = 0(r)a(t), that is, 
to solve Poisson equation for 0(r) with a fairly non-trivial right part 

V 2 <Mr,C) = -^AtGpsR^ir^PziC). (48) 
In the spherical polar coordinates we have 

j_ a r2 s^«) + ifl sin a^e) = _ 5, A> Rin ( r< _ 5/2 

r 2 Or Or r 2 sin 9 09 00 3 



Assuming a solution of the form 



(r,0) = ^P / (cos0) (50) 



and taking into account that P^(cos#) is the solution of equation 



sin 9 09 09 

we obtain 

<p 2 u(r) 5n 



1 sin6 dPt(cos9) = _ e{e+1)Pe{cos6) (51) 



r ~-^- ~ ^ + l X r ) = ~ Y Ae Gps Rl/2£/+1/2 - ( 52 ) 

The general solution of equation for (pi, which is finite at the origin, is given by 

Ur,9) = {cV + ^^T A ePs GR l /V-^P e (cos9) (53) 



Outside the star we have 



V 2 5U e = -> V 2 e (r)=O, (54) 
e = Dr- e - 1 P e (cos9). (55) 



The arbitrary constants C and D are eliminated from boundary conditions 

<Pi = 4>e 



0<pi _ 0(f) e 
r =R Or Or 



r=R 



(56) 
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which yield 



2£+V 



107T ^ 

-D = — - A e p s G 



(4£+l)(2£+l) 



R 
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Finally, we obtain 



IOtt 



3 2£ + l 



At p s G 



_ 2(2£ + 1) pl/2 £-1/2 

4£+l 



.ft (cos 0). 



(57) 



(58) 



It is worth emphasizing that following this line of argument one can get the solutions of 
Poisson equation for a more wide class of inhomogeneous star models undergoing nodeless 
spheroidal vibrations with non-rotational field of velocity. Also, the method and obtained 
solution can be useful in the study of electrodynamic problems of astrophysical interest. 



3.2. Spectral equations for frequency and lifetime of /-mode 



The mass parameter M. is given by 



R +1 p 

M = J vP (r)[al + al]dV = ^p s A 2 £ R 1/2 Jr^dr J l 2 P 2 {C) + (l - C 2 ) ( 
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(59) 



Computation of the viscous friction parameter, with the non-uniform profile of shear viscosity 
dUD, yields 

V = 2 J T)(r)aijaij dV = 2 J r)(r) (a 2 rr + a 2 ee + a 2 H + 2a 2 r9 ) dV 
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= 4vr A 2 VcR 2 "' 1 i£-l). 
The lifetime, r = 2Ai/T>, of £-pole overtone of /-mode is given by 



T f {i) = Y T »(4£ + i)(£-iy T » = -> V = J S 



(60) 



(61) 
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The integral parameter of stiffness 

1C = J [(a k V kP (r)) (a k V k U(r)) - (a k V k ) (a k V k )p(r) - p(r)(a k V fe )0(r)] dV (62) 
can be conveniently represented in the form /C = K\ + K 2 + K$, where 



J (a k V fc p(r)) (a k V k U{r))dV 
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The integral for K 2 is given by 
K 2 = - J (a fc V fc ) (a k V k )p{r) dV 
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In similar fashion, for K% we obtain 
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The resultant expression for stiffness /C reads 

?2£+1 £(4£ 2 + 2€ - 1) 



(65) 



JC = 407r2 ,2^„2 



9 



-AjGpiR 2 



(2£+l)(A£+l) 



(66) 



From analytic form of M, and /C it follows that in the singular star model under consideration 
the lowest overtone of /-mode is of dipole degree, £ = 1. 
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Fig. 1. — Frequency of fundamental modes tOf as a function of multipole degree / of global 
nodeless irrotational vibrations of singular and homogeneous star normalized to the natural 
unit of frequency ujq (upper panel), and the ratio of spectral equations for the frequency of 
fundamental mode in singular model LOf and in the Kelvin homogeneous model lok- 



The frequency spectrum of fundamental vibration mode in the singular star model reads 
2 2^(2£ + 1) — 1 Q 4tt „ GM 
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This last equation can be recast in the following equivalent form 
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(69) 



showing asymptotic behavior of the frequency at very high overtones. 



The obtained frequency spectrum of /-mode in singular star model has many features in 
common with the well-known Kelvin spectral formula for kinematically identical fundamental 
vibration mode in homogeneous star model 
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but the lowest overtone of Kelvin /-mode is of quadrupole degree, I = 2. As is demonstrated 
in Fig.l, the most essential differences between the above spectral equations are manifested 
in low-overtone domain (upper panel) and that at large £ the frequency spectra of /-mode 
in singular and homogeneous star models shear identical asymptotic behavior (lower panel). 
In the next section this last spectral formula is briefly recovered by the above expounded 
method with allow for the effect of viscous damping of /-mode whose lifetime is computed 
with non-uniform profile of shear viscosity identical in appearance to that for the profile of 
hydrostatic pressure in this model. 



4. Kelvin /-mode in the homogeneous liquid star 



In the canonical homogeneous star model of uniform density, p = constant the total 
mass has one and the same form as in above singular model, M = (4n/3)pR 3 . The internal 
and external potentials of Newtonian gravitational field are the solutions of Poisson equation 
inside and Laplace equation outside the star 
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,(71) 
(72) 
(73) 



In (172 p . W is the internal gravitational energy W of homogeneous star. The solution of 
equation of hydrostatic for pressure p(r) is given by 
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(74) 
(75) 



The last expression shows again that the star radius is determined by the central pressure 
P c related to the density p by equation of state. 

For comparison, in Fig. 2 we plot the fractional pressure and gravity potential profiles 
computed in homogeneous model 
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Fig. 2. — The fractional pressure and gravity potential profiles in singular inhomogeneous 
star model and in the canonical homogenous star model. 



and in singular inhomogeneous model 
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The governing equations for non-compressional irrotational vibrations of the homogeneous 
liquid star are 



p5v { = -Vi5p + pVi5U + V fc £7r ifc , 
Sp = -(v k V k )p, V 2 5U = 0, 

87r ik = 2r]8v ik , 5v ik = -(Vi5v k + V k 8vi), rj = r] c 



1 ^ 



The equation of energy conservation is 
d f pSv 2 



dV 



[(5v k V k )5p - p(5v k Vk)5U - 5n ik 5vik] dV. 



dt J 2 
On substituting here 

8vi(r,t) = ai(r)a(t), 6U(r,t) = (pi(r)a(t), 5p(r,t) = -(a k V k ) p(r)a(t) 



(78) 
(79) 

(80) 
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we arrive at equation of damped oscillator, M.a + Va + fCa = 0, with integral parameters 
defining the frequency and lifetime of the form 

M = J pa k a k dV, V = 2 J rj(r)a ik a ik dV, (82) 

K = - /[(aiVi)(a fc V fc )p(r)+p(a fc V fc )0(r)]dV, (83) 
Jv 

" =A?« T = p- (84) 

The only unknown quantity is the variations of the gravity potentials obeying the Laplace 
equations 

V 2 ^ = 0, V 2 5U e = (85) 
having the general solutions of the form 

SUi(r <R)= C e r e P e (()a(t), 5U e (r < R) = D t r~^ P t {Qa{t). (86) 
The arbitrary constants Cg and Dg are eliminated from the standard boundary conditions 
Ui(r{t)) + SUM*)) = U e (r(t)) + 5U e {r{t)) , (87) 

R=R 

± [UMt)) + SUMt))] = ^ Pe(r(t)) + SU e (r(t))] \ r=R (88) 

where r(t) = r[l + ct(t)Pg(cos6)]. Retaining in these equations terms of first order in a(t) 
and putting r = R we arrive at coupled algebraic equations for Cg and Dg whose solution 
leads to the following final expressions for the time-independent part of gravity potential 
(Bastrukov 1996a) 

*=2ZTll^(0> 4>e = ^GpR»*r-Wp e (0. (89) 
Computation of integral parameters of the inertia M. viscous friction V yields 

M = 4tt A 2 P R 2£+1 -4— = ;t7~7~~~T ) (90) 

V = IQn AfaR 2 ^-^-^- = 1 QnTj c R? ■ (91) 

1 ' 2i+l ' 2£+l K ' 

Following the suggestion of Cowling (1941), one can consider separately oscillations restored 
by force associated with gradient in fluctuations of pressure and force owing its origin to 
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fluctuations in the potential of gravity. In accord with this, the integral parameter of stiffness 
is written represents as a sum 

K = K p + K g , (92) 
K p = - / (a k V k ) (oiViXr) dV, K g = - [ p(a,V i )0(r) rfV (93) 

and, consequently, the frequency frequency of vibrations is also represented as a sum 



^ = ^+^, 4=** (94) 



where u p designate the frequency of p-mode and ui g is the frequency of g-mode. For K p and 
Kg we get 

16tt 2 Gp 2 R 5 16tt 2 G P 2 R 5 
p ~ 3 2£ + V 9 ~ (2£+l) 2 ■ { } 

For the frequencies of p-mode and g mode we obtain 

^ hi ^ = (96) 

v M G 9 M °2£+l v 7 

One sees that tu 2 is positive whereas u 2 is negativ^]. It follows that frequency of fundamental 
mode ojfifiQi) = [to 2 + w|] 1//2 is given by Kelvin spectral formula (1701) and for the lifetime 
of £-pole overtone, Tf(£), computed with the non-uniform radial profile of shear viscosity 
r] = rj(r) given by equation ( 1501) . we obtain 

T M = WTy T ^v- "=7 < 97 > 

It is worthwhile to compare the computed lifetime spectra with the Lamb (1945) spectral 
equation 

W) = (2f+i r ; (/ _ 1} 08) 

that has been obtained in a similar fashion but assuming that coefficient of shear viscosity 
has constant value in the entire spherical volume of homogeneous viscous liquid (see also 
Bastrukov et al 2007). It has been pointed by Jeffreys (1976), however, in the context of 



2 These sings for p-mode and g-mode are one and the same as in well-known dispersion relation of Jeans 
oj 2 = c?Ji 2 — 4nGp which characterizes propagation of sound, with velocity c s , in a self-gravitating fluid of 
constant density p (e.g. Chandrasekhar 1961, Sect. 119) 
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Fig. 3. — The fractional lifetime of £-po\e overtones of fundamental vibrational mode in 
homogeneous and inhomogeneous liquid star models with non-uniform profiles of shear vis- 
cosity in juxtaposition with Lamb spectral equation for damping time of nodeless irrotational 
oscillations of a spherical mass of viscous liquid of uniform density and constant coefficient 
of shear viscosity. 



geoseismology that the approximation of uniform viscosity does not allows for the effect of 
self-gravity on mechanical property of matter of astrophysical objects. With this in mind, 
we have supposed that it would be not inconsistent to take radial profile of shear viscosity 
similar to that for hydrostatic pressure in the star. Also noteworthy is that Newtonian law 
of shear viscosity is equally appropriate for viscous liquid and viscoelastic solid (e.g. Landau 
et al 1986). The obtained here spectral equations for the time of viscose damping of nodeless 
spheroidal vibrations may be of some interest, therefore, for the general seismology of Earth- 
like planet (Aki & Richards 2002). In Fig. 3, the obtained spectral equations for the lifetime, 
normalized to function of overtone number £ are plotted for both singular flBTl) and 

homogeneous fl97j) models in juxtaposition with the Lamb spectral formula (1981) . 



5. Summary 

The most striking differences between vibrational behavior of homogenous and inhomo- 
geneous liquid star models in the fundamental mode of global nodeless irrotational pulsations 
under the combined action of the buoyancy and gravity forces is that in the inhomogeneous 
model the oscillations are of substantially compressional character, that is, accompanied by 
fluctuations in density and the lowest overtone of /-mode is of dipole degree, as has first been 
observed in (Podgainy et al 1996). In the homogeneous star model they are characterized as 
non-compressional and the lowest overtone is of quadrupole degree. 
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Fig. 4. — Artist view of a spherical star-forming object, the stellar cocoon, whose largescale 
vibrations in /-mode can be analyzed on the basis considered highly inhomogeneous model. 

The natural unit of frequency of /-mode, ujq = a/ GM/R 3 , has one and the same 
form as that for gravity or g-modes. These letter modes are of particular importance for the 
asteroseismology of transitory pre-white dwarfs objects and young white dwarfs (e.g. Hansen 
& Kawaler 1994) whose variability of electromagnetic emission is attributed to non-radial 
gravity- driven seismic vibrations (Koester & Chanmugam 1990, Fontaine & Brassard 2008, 
Winget & Kepler 2008). To this end, it is worth emphasizing, the existing modal classification 
of the g-mode spectra presumes that white- dwarf-forming object oscillates in the standing- 
wave regime in which the frequency is determined by nodal structure of fluctuating variables 
(e.g. Unno et al 1989, Winget & Kepler 2008). According to standard nomenclature, the 
g-modes of standing-wave regime of non-radial pulsations are specified as n gt, where n > 1 is 
the node number and i is the multipole degree of overtone. In the meantime, the regime of 
nodeless oscillations, designated as oge, remains less studied. Also, the considered model can 
be invoked in assessing variability of emission from star-forming clouds like stellar cocoon, 
pictured in Fig. 4, as caused by its vibrations in fundamental mode. Motivated by these 
arguments we have investigated an analytic model of a highly inhomogeneous star undergoing 
global nodeless irrotational vibrations. In work (Bastrukov 1996b) these kind of nodeless 
oscillations have been studied in the model of homogeneous liquid layer. In recent papers 
(Bastrukov et al 2007b, 2008), the nodeless vibrations of a solid star driven by restoring 
force of elastic deformations and locked in the peripheral finite-depth mantle, has been been 
studied in the context of seismic vibrations of the neutron star crust in the approximation 



-18- 



of uniform density of crustal matter. Understandably that in real white dwarfs and neutron 
stars the bulk density is non-homogeneous function of position. However, this realistic case 
of nodeless irrotational vibrations entrapped in the peripheral finite-depth seismogenic layer 
with non-uniform density profile demands more elaborate mathematical treatment which 
will be the subject of our forthcoming paper. 

The authors are grateful to Gwan-Ting Chen (NTHU, Hsinchu, Taiwan) for helpful 
assistance. This work is a part of projects on investigation of variability of high-energy 
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